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ABSTRACT 

The Frying Pan (G315.9— 0.0) is a radio supernova remnant with a pecuhar hnear feature 
(G315.78— 0.23) extending 10' radiahy outward from the rim of the sheh. We present radio imag- 
ing and polarization observations obtained from the Molonglo Observatory Synthesis Telescope and 
the Australia Telescope Compact Array, confirming G3 15. 78— 0.23 as a bow-shock pulsar wind nebula 
(PWN) powered by the young pulsar J1437— 5959. This is one of the longest pulsar tails observed in 
radio and it has a physical extent over 20 pc. We found a bow-shock stand-off distance of 0.002 pc, 
smallest among similar systems, suggesting a large pulsar velocity over 1000 km s~^ and a high Mach 
number ~200. The magnetic field geometry inferred from radio polarimetry shows a good alignment 
with the tail orientation, which could be a result of high flow speed. There are also hints that the 
postshock wind has a low magnetization and is dominated by electrons and positrons in energy. This 
study shows that PWNe can offer a powerful probe of their local environment, particularly for the 
case of a bow shock where the parent supernova shell is also detected. 

Subject headings: ISM: individual (G315.78— 0.23, G315.9— 0.0) — ISM: supernova remnants — pul- 
sars: individual (PSR J1437— 5959) — radio continuum: ISM — stars: neutron — 
stars: winds, outflows 



1. INTRODUCTION 

Neutron stars are born with large space velocities of 
a few hundred kilometers per second. They eventu- 
ally escape their natal supernova remnants (SNRs) af- 
ter ^ 10^ yr and travel supersonically through the in- 
terstellar medium (ISM), which typically has a sound 
speed of 1-10 km s~^ in the cold and warm phases. In 
this case, the pulsars' relativist ic outflows are confined 
by ram pressure of the ISM, resulting in bow-shock 
pulsar wind nebulae (PWNe). These systems exhibit 
comet ary morphologies with long tails, and emit broad- 
band synchrotron radia tion from radio to X-rays (see, 
iGaensler fc Slanell2006L for a review). As compared to 
PWNe within SNRs, bow-shocks are governed by sim- 
pler boundary conditions, therefore, providing a powerful 
probe of the local ISM environment. 

At radio frequencies, PWNe generally show a larger 
extent than their X-ray counterparts, due to long syn- 
chrotron cooling timescale of radio-emitting particles, 
which can be up to 10^-10^ yr. Therefore, radio PWNe 
can act as direct calorimeters to reflect the integrated 
history of the systems. Table [1] lists the longest pul- 
sar tails detected in radio. The best example is the 
Mouse PWN, its morphology consists of a broad cone- 
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shaped outer region and a collimated narrow inner re- 
gion, which were suggested to be postshock flows from 
the for ward and backward t ermination shocks, respec- 
tively (jGaensler et al.l l2QQ4f ). Radio observations of 
PWNe also reveal the nebular magnetic field structure 
via polarimetry. Previous studies indicate a diverse B- 
field configuration among bow-sh ock PWNe, from a heli- 
cal field perpendicular to the tail (Ng et al."2010') to a B- 
field parallel to the tail (Yusef-Zadeh & Gaensler 20051), 
but the origin of the variation remains a puzzle. 

The "Frying Pan" SNR (G315.9-0.0) is a faint 15'- 
diameter radio shell with an intriguing 10' linear protru- 
sion, 0315.78— 0.23 (the "handle"), extending northwest 
from the rim, such that the overall morphology resembles 
a Frying Pan (see Figured]). Based on the lack of an IR 
counterpart, W^hiteo ak fc Green (1996) concluded that 
the radio emission is nonthermal and the source is an 
SNR. In this case, we suggested that the highly unusual 
linear feature is due to the trailing relativistic wind of a 
fast pulsar moving outward from the central birth place. 
We have recently performed a deep radio search near the 
tip of the protrusion and discovered a new radio pulsar, 
J 143 7— 5959, providing s trong support to the above pic- 
ture (|Camilo et al.|[2QQ9l ). This pulsar has a short period 
of 61.7ms, a characteristic age of 110 kyr, and a high 
spin-down power ^=1.5 x lO^^ergs"^. While its large 
dispersio n measure of 549pc c m~^ suggests a distance of 
8kpc, as iCamilo et al.l (|2009[ ) noted, the estimate is un- 
certain and could be easily off by 25%. 

In this paper, we focus on the handle of the Frying Pan 
and establish its PWN nature using new and archival 
radio observations from the Molonglo Observatory Syn- 
thesis Telescope (MOST) and the Australia Telescope 
Compact Array (ATCA). Its projected length over 20 pc 
makes it one of the longest pulsar tails observed. We 
show that this unique object can provide crucial infor- 
mation on the pulsar environment that is not easily ac- 
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TABLE 1 

Long Pulsar Tails in Radio, Sorted in Physical Length 



PWN 


SNR 


PSR/Neutron Star 


Length 


Velocity 


Distance-*^ 


Reference 








(pc) 


(kms-i) 


(kpc) 




G315.78-0.23 


Frying Pan 


J1437-5959 


20 


-1000 


8 


This work 


Mouse 


J1747-2958 


16 


300 


5 


1,2 


G319.9-0.7 




J1509-5850 


10 


300 


3 


3,4 


G68.77+2.82 


CTB 80 


B1951+32 


5 


274 


2 


5,6 


G34.56-0.50 


W44 


B1853+01 


2 


375 


2.6 


7 


Duck 




B1757-24 


1.5 


<360 


5.2 


6 


G47.38-3.88 




B1929+10 


1.3 


177 


0.36 


8,9 


G189.23+2.90 


IC 443 


J061705. 3+222127 


1 


230 


1.5 


10,11 


N157B2 


N157B 


J0537-6910 


30 


? 


50 


12 


G309.92-2.5l3 




J1357-6429 


7 


? 


2.5 


13 



Ref eren ces. — (plGaensler et aLl (120041): (2) [H ales et al. (2009); (3) Kargalts ev et al.l 
(I2OO8II: (4 ) iNs et aLl (12010; (5) 'Castelletti et al. (2003); (6) Zeieer et al. (200^); (7) 



iFrail et aLl (1996); (8) Chatteriee et al. (2004); (9) Becker et al. (2006); (10) Olbert e tall 
(|200ll ): (11) ..Gaensler et al. (2006); (12) Lazendic et al. (2000 ): (13) Chang et aL (2011). 

^ All distance estimates are subject to large uncertainties, except for PSR B1929+10, which 
is obtained from parallax. 

^ This syst em is more likely to be a PWN crushed by the SNR reverse shock, instead of a 
bow- shock ( Ivan der Swaluwl[2003 : IChen ^Tall [20061 ). 
3 The radio counterpart is not confirmed. 



cessible in ottier systems. 

2. OBSERVATIONS AND DATA REDUCTION 

We carried out new radio imaging observations of the 
Frying Pan with the ATC A at 3 and 6 cm in four differ- 
ent array configurations, and processed archival ATC A 
data at 13 and 20 cm, incl uding those from the S outh- 
ern Galactic Plane Survey (jHaverkorn et al.ll2QQ6l ). Ta- 
ble [2] lists the parameters of all datasets used in this 
study. The last two observations in the list were car- 
ried out after the Compact Array Broadband Backend 
(CABB) upgrade (jWilson et al.l 1201111 . which provides 
2 GHz bandwidth, nearly 20 times improvement over the 
pre-CABB era. Note that with such a large bandwidth, 
the 3 and 6 cm data and images we refer to throughout 
this paper actually cover 3.0-3.7 cm and 4.7-6.6 cm, re- 
spectively. The pre-CABB observations at 3 and 6 cm 
were made with a two-pointing mosaic to cover 5' from 
the tip of G315.9— 0.0. One more pointing was added in 
the CABB observations to extend the coverage over the 
entire 8'-handle. The fiux density scale is set by obser- 
vations of the primary calibrator, PKS B1934— 638. A 
secondary calibrator, PKS B1414— 59, which is 2?4 away 
from our source, was observed every 30 minutes to deter- 
mine the antenna gains. 

All data reduction was performed using the MIRIAD 
packag43- After discarding the edge channels and chan- 
nels affected by self-interference, we obtained a usable 
bandwidth of 104 MHz and 1848 MHz for the pre-CABB 
and CABB data, respectively, at each wavelength. The 
longest baselines in the 13 and 20 cm data were excluded 
to provide a uniform u-v coverage. We examined the 
visibility data to reject bad data points due to poor at- 
mospheric phase stability, then applied gain, bandpass, 
fiux, and polarization calibrations. Mosaiced intensity 
maps were formed using the multifrequency synthesis 
technique with uniform weighting at 13 and 20 cm and 
with natural weighting at 3 and 6 cm. The maps are de- 

^ |http : //www . atnf . csiro . au/comput ing/sof tware/miriad/ 1 



convolved using a maximum entropy algorithm (MOSMEM; 
iSault et anil999f ) and restored with Gaussian beams of 
FWHM 23'' X 20'' at 20 cm, 12" x 9'.'7 at 13 cm, and 
3'.'7 X 3'.'3 at both 3 and 6 cm. The final maps have rms 
noise of 1, 0.1, 0.011, and 0.016 mJybeam-^ at 20, 13, 6, 
and 3 cm, respectively. While the last two values are very 
close to the theoretical limits, the 20 cm map is severely 
contaminated by a strong source nearby and the noise in 
the 13 cm map is also slightly higher than the theoretical 
level of 0.06 mJy beam~^. 

For polarimetry, we only focused on the 3 and 6 cm 
data since the other observations lack the sensitivity and 
suffer from beam depolarization. To mitigate the prob- 
lem of bandwidth depolarization due to the large frac- 
tional bandwidth of CABB data, we divided the 6 cm 
CABB data into 33 sub-bands (56 MHz each) and the 
pre-CABB data into 2 sub-bands. The problem is less 
severe at 3 cm, therefore, we only divided the CABB data 
into 5 sub-bands. Individual Stokes Q and U maps for the 
sub-bands are formed using the same procedure as above 
for Stokes I. We attempted t o apply the rotation measure 
(RM) synthesis technique (Brentiens fc de Bruvnl [20051 ) 
to determine the foreground Faraday rotation. Although 
the result i s generally co n sisten t with the pulsar RM re- 
ported by iCamilo et al.l (|200 9l). the uncertainty is too 
large to be useful. Therefore, we adopted the pulsar 
value of —700 ± 25radm~^ to derotate the individual 
Stokes Q and U maps. The error in RM corresponds to 
uncertainties of 1?5 and 4? 3 in the polarization angles at 
3 and 6 cm, respectively. Finally, the derotated maps in 
each frequency channel were weighted according to their 
rms noise and co-added to generate maps of polarized 
intensity and position angle. 

In addition to the radio observations, we also pro- 
cessed the Fermi Large Area Telescope data to search 
for pulsations in 7-rays. We carried out the analysis us- 
ing the Fermi Science Tools v9rl8p6 and considered only 
the "diffuse" class of events with a cut at zenith angles 
>105°. We only used the first year of Fermi data taken 
before 2009 September 27, when the pulsar ephemeris 
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TABLE 2 

Parameters for the ATCA Observations of G3 15. 9 -0.0 



Obs. Date 


Array 


Center Frequency 


Usable Bandwidth 


No. of 


Integration 




Configuration 


(MHz) 


(MHz)i 


Channels -"^ 


Time (hr) 


1998 Jul 5 


750E 


1384, 2496 


104 


13 


11 


1998 Dec - 2000 Aug^ 


various-^ 


1384 


104 


13 


0.33 


1999 Mar 25 


1.5B 


1384, 2496 


104 


13 


12 


1999 Apr 20 


1.5C 


1384, 2496 


104 


13 


12 


2008 Jul 23 


H214 


1344, 1436 


104 


13 


9 


2008 Dec 8 


750B 


4800, 8640 


104 


13 


11.5 


2009 Feb 18 


EW352 


4800, 8640 


104 


13 


12 


2009 Jun 13 


6A 


5500, 9000 


1848 


1848 


12 


2009 Jul 31 


1.5A 


5500, 9000 


1848 


1848 


12.5 



^ Per center frequency. 

2 Part of the Southern Galactic Plane Survey (|Haverkorn et al.ll2006l l. 



given bv lCamilo et al.l (|2009D is valid. We extracted 0.1- 
lOGeV photons within 1° of the pulsar position, ap- 
plied the bary centric correction to their arrival times, 
and folded them according to the radio ephemeris. How- 
ever, no pulsation is detected. This is not unexpected 
given the large pulsar distance. 

3. RESULTS 
3.1. Overall Morphology and Flux Densities 

Figure [1] presents radio continuum images of the Fry- 
ing Pan at 36, 20, 13, 6 and 3 cm. All of them were 
taken with the ATCA as described in Section [2] ex- 
cept the 36 cm one, whi ch is from the MOST Sup ernova 
Remnant Catalogu^B (jWhiteoak fc GreenI I1996D . The 
13 cm image has be en reported in our previous paper 
(|Camilo et al.ll20Q9l ). The radio shell exhibits a circu- 
lar morphology with inner and outer radii of 4. '5 and 7', 
respectively. The overall shell is brighter to the west, 
particularly in the northwest where it intersects with the 
pulsar tail, but it is very faint and diffused to the north- 
east. Excluding the handle, the shell has flux densities 
of 0.4 ± 0.1 and 0.5 ± O.lJy at 36 and 13cm, respec- 
tively. The image at 20 cm is severely contaminated by 
strong sources nearby, precluding any detailed flux mea- 
surement. Also, the 3 and 6 cm observations were de- 
signed to map the tip of the linear protrusion, hence, 
they are insensitive to the overall SNR structure. 

As shown in Figure [TJ the handle is highly collimated 
and extends 8. '5 from the shell to the pulsar, at a posi- 
tion angle of 314° (north through east). It has a width 
of 6'' at the tip, which is only resolved in the 3 and 6 cm 
images, then gradually widens to ~ 60'' at the intersec- 
tion with the shell. At the source distance of 8 ± 2 kpc, 
these correspond to a physical length of 20 ± 5 pc and 
width from 0.23 ± 0.06 pc to 2.3=b0.6pc. Beyond 5' from 
the pulsar, the 3 cm emission falls below the observation 
sensitivity. We note that extrapolating the handle inside 
the remnant slightly misses the shell's geometric center 
by ~1.'5 to the northeast. 

The entire handle has flux densities of 110 ±40, 35 ± 1, 
26 ± 10 and 18 ± 5mJy at 36, 13, 6, and 3 cm, respec- 
tively. These measurements have large uncertainties due 
to strong variation in the diffused background. There- 
fore, we did not attempt to derive the radio spectra index 
from them, instead, we employed the spectral tomogra- 

|http : //www . physics . usyd . edu . au/sif a/Main/MSC | 



phy technique between the 3 and 6 c m data. We followed 
a simi lar procedure as described by iGaensler fc Walla"c3 
(|2003f ) and smoothed the maps to 10'' resolution to boost 
the signal-to- noise ratio. We found a spectral index 
a = — 0.45 ± 0.1 (defined as oc i/*^), but no signifi- 
cant spatial variation of a. 

3.2. The Handle 

Figures [2ja) and (b) show the 3 and 6 cm intensity 
maps zoomed-in near the pulsar. The detailed morphol- 
ogy of the handle is very similar at both wavelengths. 
There is a hint of unresolved compact emission at the 
pulsar position at 6 cm (see the contours in Figure [2fc)) 
with a flux density 60 ± 30/iJy. However, this is not ob- 
served at 3 cm. The radio emission of the handle peaks 
at 45" from the pulsar, with flux densities of 0.25 and 
0.20 mJy beam" ^ at 6 and 3 cm, respectively. Further 
away, there is a kink at 1' from the pulsar, such that 
the handle changes orientation by 34° and runs towards 
east for 10", then switches back to the original direction. 
The reverse occurs around 30" further southeast. As a 
net result, a 40"-long section is displaced by 6" to the 
northeast. 

For polarimetry results, the polarization 5- vectors, af- 
ter correcting for foreground Faraday rotation, are over- 
plotted in Figures [2fa) and (b). These indicate a highly 
ordered magnetic field structure. The 5-field well aligns 
with the handle's elongation, and changes direction along 
with the handle across the kink. Farther from the kink, 
the field orientation remains constant until the polar- 
ization signal falls below detection beyond 2. '4 from the 
pulsar. Figures [2jc) and (d) show the polarized intensity 
maps at 6 and 3 cm. The polarized emission generally 
follows the total intensity, except across the kink, where 
it shows a slight depolarization. We estimate that the 
polarization fraction is up to ~40% and ~50% at 6 and 
3 cm, respectively. A reduced polarization fraction at 
lower frequenc y has been observed in other sources (e.g., 
iNg et alll2010f ) and it could be attributed to beam depo- 
larization. However, we note that these values are rough 
estimates due to high background. Also, variations in po- 
larization angle could shift power in the polarized image 
into smaller scales, to which the interfero metric observa- 
tions are more sensitive (see discussion in IGaensler et al.] 
Il999f ). Therefore, the absolute polarized flux could be 
biased and the maps should be taken as relative mea- 
surements only. 
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Fig. 1. — Radio continuum images of G315.9— 0.0 at different wavelengths. The 36 cm image is adopted from the MOST Supernova 
Remnant Catalogue ( Whiteoak & Green 1996) and the rest are taken with ATCA. The large and small boxes in the 36 cm image indicate the 
field of views of the 6 and 3cm images, respectively. Position of PSR J1437-5959 [(J2000.0) R.A.=14h37™01?91(3), decl. = -59°59'01^'4(3) 
(jCamilo et al.|[2009l )] is marked by the cross. The restoring beams are shown at the lower right of each panel. The gray scales are linear 
with the scale bars in units of Jybeam"-*^. The 13 cm map is same as the one reported in' ICamilo et al.j ((20091) . 
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Fig. 2. — Top: radio continuum images of the tip of G315.9— 0.0 at: (a) 6 cm and (b) 3 cm, overlaid with the corresponding polarization 
5- vectors that show the intrinsic magnetic field orientation. The gray scales are linear, ranging from —0.045 to 0.22 mJy beam" -"^ at 6 cm 
and from —0.02 to 0.3 mJy beam" -"^ at 3 cm. The vector lengths are proportional to the polarized intensity, with the scale bars at the lower 
right representing 0.1 mJy beam" Bottom: linear polarized emission at: (a) 6 cm and (b) 3 cm, overlaid with total intensity contours at 
levels of 0.03, 0.073, 0.15, and 0.22 mJy beam"! for 6cm and 0.06, 0.12, and O.lTmJy beam"! for 3 cm. The gray scales are linear, ranging 
from to 0.15 mJy beam" -"^ at 6 cm and from to 0.2 mJy beam" -"^ at 3 cm. In all panels, position of PSR J1437— 5959 is marked by the 
cross and the restoring beams are shown at the lower right. 



4. DISCUSSION 

4.1. Nature of the Handle and Its Association with the 

SNR 

Our results show that the handle of the Frying Pan 
(G315.78— 0.23) has a high degree of linear polarization 
with a flat radio spectrum, hence confirming that it is a 
PWN powered by PSR J1437-5959. If the source is at 
a distance of 8 kpc, then it would have a physical length 
over 20 pc, the longest pulsar tail ever detected. The 
source morphology suggests that it is a bow shock with 
the pulsar traveling highly supersonically in the ISM. 
As radio PWNe generally have a very long lifetime (see 
Equation (p!Q|) below), the fact that the tail points toward 
the general direction of the remnant center argues for an 
associati on between PSR J1437-5959 and G315.9-0.0 
(but see. lZeiger et al.|[2QQ8[ ). We performed Monte Carlo 



simulations to quantify the probability of chance asso- 
ciation. A^^_eni2la^^ed_j^j;ea^ synthesis 
model (|Faucher-Giguere" fc KasDil[2QQ6l : A. Bouchard et 
al. 2012, in preparation) and found that at this sky po- 
sition, an average of 0.2 field radio pulsars are expected 
within 20' from the shell, younger than 10^ yr and with 
a radio flux larger than that of PSR J1437-5959. If 
we further require the pulsar's projected motion in the 
past to intercept any part of the shell, then the number 
greatly reduces to 0.0012, therefore, strongly suggesting 
an association between PSR J1437— 5959 and the Frying 
Pan. Although the tail does not point toward the exact 
geometry center of the shell, this is not uncommon and 
could be the result of asymmetric supernova explosion, 
the progenitor's prop er motion, or a sligh^_gradient in 
the ISM density (se e, iGvaramadzg 120061 : iNg et al.l 120071 : 
IVigelius et"alll2007[ ). 
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To estimate the kinematics of the system, we extrap- 
olate the tail inside the shell and take the closest ap- 
proach to the center as the supernova site, which is 14.'5 
away from the pulsar location. This implies a distance 
/ = 1.0 X 10^^ (ig cm traveled by the pulsar since birth, 
at the source distance of Sdgkpc Here we assume the 
pulsar motion lies in the plane of the sky for simplic- 
ity. The pulsar velocity and age t are then related by 
= l/t. Even if the pulsar travels as fast as 2000 km s~^, 
the system would still be older than 15kyr, such that 
G315.9— 0.0 has long passed the free expansion phase. 
On the other hand, the transition from the Sedov phase 

to the snowplow phase occurs at ttr ^ 29E^(^^nQ^^^^ kyr 

5/17 T / 

with a shell radius i?tr ~ 19.1£^5{ Uq , where E^i is 
the supernova explosion energy in units of 10^^ erg and 
no is the ISM number density in cm~^ (|Blondin et al.l 
|1998[ ) . Our observed shell radius of 7' gives a physical size 
R = 16.3 (igpc, slightly smaller than i?tr for the typical 
values of £^51 = 1 and no between 0.1 and 1. Therefore, 
we adopt the Sedov solution 



R= 1.15 



P 



1/5 



(1) 



where Eq 
gcm~^. 



10^^£^5i erg and p is the ISM density in 



4.2. 



Pulsar Velocity and Age 

The bow-shock stand-off distance rbs ahead of the 
PWN is determined by the ram pressure balance between 
the pulsar wind and the ISM. For isotropic wind. 



(2) 



Assuming a constant p across the field, we can solve 
Equations ([1]) and (|2]) to obtain 



E 



1/2 



AncEoP 



R 

LIS 



5/2 



0.0024 dg^^^sYpc . 



1/2 

This corresponds to an ang ular scale O'.'Oe^^ K 



51 



(3) 
too 

small to be resolved by the observations. It is interest- 
ing to note that this value is independent of the sys- 
tem age, pulsar velocity and ambient density. Also, it 
is one of the smallest stand-off distances among bow- 
shock PWNe, e.g., an order of magnitu de smaller than 
that of the Mouse (|Gaensler et al.ll2QQ^ . This indicates 
either a very large Mach number or a very dense ambi- 
ent medium. We argue that the latter case is less likely, 
since the CO map (jPame et al.ll2001[ ) shows no evidence 
of molecular clouds in the region. As mentioned, there 
is a marginal detection of a compact feature at the tip 
of the PWN with a flux density 60 zb 30/iJy at 6 cm. If 
this is from the pulsar, then c omparing; to the p ulsed 
flux density of 75/iJy at 20 cm ( Camilo et al.l 12009) in- 
dicate a rather flat radio spectrum. As an alternative, 
this could also be the surface of the backward termi- 
nation shock, which has a n extent of 5rhs acco rding to 
hydrodynamic simulations (jGaensler et al.ll2004f ). such a 
small-scale structure would not be resolved in our im- 
ages. Further observations with high time resolution 



could distinguish between these cases. Finally, we note 
that the PWN half- width near the tip is about 40rbs, 
much greater than the val ue found in magnet ohyd rody- 
namic (MHD) simulations (jBucciantini et ani2005f ). The 
discrepancy could possibly be due to anisotropy of the 
pulsar wind or mass loading through co ntaminati on (se e 
IVigelius et alll2007l : iRomani et al]l2010l : lLvutikovll2003[ ) . 

Assuming cosmic abundances, we substitute rbs into 
Equation ^ to obtain 



= 1700 Eli^d^ 



-3/2 -1/2 
^0 



kms 



and hence 



t = 19K 



-1/2 .5/2^1/2 



51 



kyr . 



(4) 
(5) 

While we argue above that no cannot be too high, no 
cannot be too low either because of the detection of the 
radio shell in G315.9— 0.0, hence, we can rule out the hot 
ISM phase. For a typical value of no between 0.1 and 
1 in the warm ISM phase, the inferred velocity would 
be nearly 2000kms~^with a young system age in the 
order of 10 kyr, as proposed bv I Camilo et al.l (|2009l ) . This 
indicates t < ttr, providing support to the Sedov phase 
for the SNR. We also note that the age estimate is much 
younger than the pulsar's characteristic age of 110 kyr, 
which could be reconciled if the pulsar was bo rn with a 
spin period close to the present-day value (see iNg et al.l 
120071 ). With these results, we can estimate the pressure 
^head at the head of the bow shock and the pulsar Mach 
number A4 using 



pvl = 6.8x10"^ £;5i^^8^<^yncm"^ 



(6) 

where Pism is the ambient pressure and 7 = 5/3 is the 
ISM adiabatic index. Since different phases of the ISM 
are generally in pressure equilibrium, Pism should be 
more uniform than p, hence, the Mach number estimate 
would be more certain than the pulsar velocit y estimate 
(see the detailed discussion by Kargalt sev et al.ll2008f ). 
A typical ISM pressure of 10~^^dyncm~^ (correspond- 
ing to Pisu/k = 7250cm-3K) gives M = 200 E^'^d'^^'^ , 
suggesting a highly superso nic pulsar niotion. Fo r 
smaller values of Pism (e.g., iFerrierd [20QTI : ICo^ 120051 ). 
the inferred Mach number would be even higher. While 
such a highly supersonic motion is expected to produce 
a very sharp Mach cone with a small half-opening an- 
gle ~ sin~^ (1/200), the observed width of the tail in- 
creases from 6" to 60^' over its entire length of 7\ im- 
plying ^ sin~^(l/15). The discrepancy could be the 
result of overpressure in the tail, such that 

Ptaii ^ p{v* sin 6>)^ = Phead siu^ 

= Phead/15'= 0.004Phead • (7) 

See Section |431 below for further discussion. 

4.3. The Kink 

Figure [2] shows a kink in the tail at V from the pul- 
sar. We argue that this is unlikely caused by flow in- 
stability or by a global ISM pressure gradient, since the 
magnetic field is highly ordered across the kink and the 
tail resumes its original orientation further downstream. 
Instead, it could be the result of interstellar turbulence. 
Comparing the kink's 6''-lateral displacement to its V- 
separation from the pulsar , we can deduce a turbulent 
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velocity Vt = {&^/V)v^ = 0.1'^;*, which is much higher 
than the typical sound speed of ~10kms~^ in the warm 
ISM phase, giving evidence of supersonic turbulence. As 
this PW N system has a rather high flow speed (see Sec- 
tion |43] below) , the postshock wind could become com- 
pressible. The compression by the ISM turbulence would 
result in enhanced synchrotron emissivity, which could 
naturally explain the coincidence of the kink and the ra- 
dio peak. Also, it would introduce disorder in the mag- 
netic field, reducing the polarized fraction near the kink 
as observed. 

4.4. Magnetic Field Strength 

To estimate the magnetic field strength in the PWN, 
we define km = Us/Up to be the magnetization of the 
postshock wind, i.e., the ratio between the magnetic 
and relativistic particle (electron+ion) energy densities 
(see lPavlov et al.ll2003l ). MHD simulations show a rela- 
tively uniform mag:n etization in the tail wi th km of the 
typical order of 0.1 (Buccian tini et aP 1200 5). Let rj be 
the ion to electron energy density ratio, we then have 
Up = {1 -\- r])Ue and we can estimate the average mag- 
netic field strength B from standard synchrotron theory 



Ub 



fcm(l +7?) 



Cl2 



(8) 



where ci2 is a constant that depends weakly on the spec- 
tral index and on the lower and upp er frequency limits of 
the emission (see iPacholczvkl [19701 ) and L/V is the syn- 
chrotron volume emissivity. Assuming a simple power- 
law spectrum between 10^ and 10^^ Hz as in Ng et al. 
(|2010l ). our measured fiux density and spectral index 
yield 



B 



87rci2km{l^v) ^ 



2/7 



(9) 

We note that this result is insensitive to the frequency 
limits, and to km and 77. For example, changing the upper 
limit to 10^^ Hz only lowers B by 13%. If we Ignore 
any ions (i.e., r] = 0) and take km = 0.1, we obtain 
B = 30/iG. Adopting km = 0.01 yields B = 16/iG, while 
the equipartition case {km ^ 1) gives B = 60/iG. For 
emission at 6 cm {u ^ 5 GHz), the synchrotron lifetime 
is 



ISB 



-3/2 -1/2 



kyr : 



540 k-^^\l^r])-^^^ dl^^kyi , 



(10) 

sufficiently longer than the estimated age of the sys- 
tem (see Equation (|5|)). Comparing to the synchrotron 
cooling, inverse Compton scattering with the cosmic mi- 
crowave background is negligible. Therefore, the radio 
PWN should trace the pulsar motion over a long time, 
providing support to the association with G315.9— 0.0. 

4.5. Pressure and Flow speed 

The magnetic field estimate can give us a handle on 
^ISM, which consists of the magnetic pressure and the 
particle pressure. Taking the pulsar wind as relativistic 
ideal gas, the particle pressure is given by Up/ 3, hence, 

Ptau=C + ^ = (i + 3i;:)C = i-5xio-- 

{l + ^)kT{l + v)^^'ds''^'dyncm-^ . (11) 



This confirms that the system is highly overpressured as 
compared to the environment. To match Ptaii we derived 
in Equation ([7]), it would require 77 < 1 and km ^ 0.02, 
hence, this gives some hints that the postshock wind is 
dominated by electrons and positrons with a small mag- 
netization. Finally, we note that Ptaii above is much 
lower than the value 0.02Phea d found in MHD simula- 
tions (|Bucciantini et aP l2005f ) , although those simula- 
tions only extend for a few rbs in the tail direction. 

We can estimate the average fiow speed Vf of the post- 
shock wind based on energy conservation. If we assume 
the pulsar spin-down energy all converts into the mag- 
netic and particle kinetic energy downstream, then 

E = A{vf + v^){Ub + Up) = A{vf + V*) ( 1 + 

thm f 

(12) 

where A is the tail's cross-sectional area. We modeled 
the tail with a truncated cone of diameter from ^" to 
60'', then A scales with the distance z (in arcminutes) 
from the pulsar as A{z) « 4.0 x 10^^(1 + l.lzf'^ cm~^, 
giving 



Vf ^ 2.3x10^(1+1.1^) 



^(l+r?)-^/^d3-^°/^kms-\ 

(13) 

here is ignored as it is negligible. Our result indicates a 
rather high flow speed near the pulsar {z < 1): for r] = 
and /c^=0.1, we obtain Vf ~ 0.06c. If A;^=0.01, then Vf 
is a factor of 2.5 lower. These numbers are slightly higher 
than those observed in other bow-shock PWN sy stems 
fe.g.. iGaensler et al.ll200l iKargaltsev et al.l 120081), but 
lower than the simulation prediction (jBucciantini et aP 
I2OO5I ). We note that while synchrotron loss is negligible, 
other energy dissipation mechanisms, e.g., induced by 
flow instability, mass- loading or ions in the wind, coul d 
signiflcantly reduce the flow speed (see iLvutikovl |2QQ3[ ) . 
Combining our radio results with future X-ray observa- 
tions can provide a more complete picture of the particle 
spectrum, thus, better accounting for the total energet- 
ics of the wind to reveal the flow structure. Spatially- 
resolved X-ray spectroscopy can also indicate any de- 
celeration in the flow. For particles emitting in the X- 
ray range {u ^ 5 x 10^^ Hz), Equation (p!Q|) gives a syn- 
chrotron cooling time of 150 yr. Hence, a flow speed of 
0.06c would result in a 3 p c-long X-ray PWN, com parable 
to other long X-ray tails (jKargaltsev et al.ll2008f ). 

4.6. Comparison with Other Long Pulsar Tails in Radio 

In Table [H we summarize the physical extents and 
space velocities of the longest pulsar tails detected in 
radio. We found only a very week correlation between 
these two quantities. This is not surprisingly, since we 
expect that the radio emission should also depend on 
the physical conditions of the local environment, such as 
ambient density and magnetic fleld. The three longest 
radio tails: the Frying Pan handle (20 pc), the Mouse 
(16 pc) and G319.9-0.7 (10 pc), ah ha ve very high Mach 
numbers (200, 60 and 30, respectivelv: iHales et al.ll2009l : 
IKargaltsev et al.ll2008f ) and exhibit smooth morphologies 
with highly ordered magnetic flelds. This indicates a 
smooth flow structure in these systems, which helps ex- 
plain the high degree of collimation over long distances. 
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On the other hand, these three sources show some dis- 
crepancies in the detailed morphology and magnetic field 
structure. The radio emission of the Mouse peaks at 
the head then gradually fades downstream, and it shows 
a cone-shaped outer region and a narrow inner region. 
For the handle of the Frying Pan, the tip is very faint 
and a two-component structure is not obvious, although 
the unresolved source at the pulsar position may corre- 
spond to the inner region, or the outer region could be 
too faint to detect. Both the Frying Pan handle and 
the Mouse have a very similar magnetic field geometry 
that runs parallel to the tail (Yusef-Zadeh & Ga enslerl 
I2QQ5D . whil e G319.9-Q.7 s hows a 5-field perpendicular 
to the tail (jNg et al.ll2QTQl ). The origin of this diversity 
is unclear. It may depend on the pulsar spin orientation 
relative to its motion or the fiow condition of the system. 
In the regime of high postshock fiow speed, the magnetic 
field may be sheared by the fiow so that the field lines 
are stretched along the tail. More examples are needed 
to verify this picture. Finally, we note that the polar- 
ization vectors in both the Mouse and G3 19. 9— 0.7 ex- 
hibit abrupt switches at large distances from the pulsar 
^Yusef-Zadeh & Bally 1987; Ng et al. 2010), which is not 
observed here. Therefore, we believe that this is not a 
universal feature among long pulsar tail. 

5. CONCLUSIONS 

In this paper, we have presented a detailed radio 
study focusing on the handle of the Frying Pan SNR 
(G315.9-0.0) using ATCA and MOST observations. Our 
results confirm the 20 pc-long protrusion as a bow-shock 
PWN associated with PSR J1437-5959. The detection 
of both the PWN and the parent SNR makes this sys- 
tem a powerful probe of the pulsar environment. The 
very small bow-shock stand off distance of 0.002 pc we 
deduced implies a large pulsar velocity over 1000 km s~^ 
and a high Mach number ~200. The PWN has a highly 
ordered magnetic field that runs parallel to the tail ori- 
entation and has an average strength ~20/iG. Based on 
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the pressure estimate and the tail's opening angle, we 
found some hints of a small magnetization in the post- 
shock wind of the order of 0.01-0.1. Also, the ions in 
the wind could be less energetic than the electrons and 
positrons. Finally, we suggest that the kink observed in 
the tail is a result of supersonic ISM turbulence. 

As a caveat, the main uncertainties in our derived 
quantities come from the distance estimate. It is not 
uncommon for dispersion measure-based distances to be 
off by 25% or more. While we may have to wait for the 
Square Kilometer Array for a direct distance measure- 
ment with parallax (Smits et al.ll2011f ). multiwavelength 
studies are essential to further understand the physics of 
this remarkable system. In particular. X-rays observa- 
tions can provide a complementary picture of the parti- 
cle energetics and reveal the flow structure of the high- 
est energy particles. If TeV emission is detected, as in 
other PWNe, a joint modeling of the synchrotron and in- 
verse Compton scattering emissions can give an accurate 
measurement of the magnetic field strength. For a pul- 
sar space velocity of 2000 km s~^, a direct proper motion 
measurement with the ATCA would require a time span 
over 10 yr. Therefore, other techniques, such as vey long 
baseline inter fer mo metry or timing observations, would 
be more practical. The latter is also essential for a deeper 
search of 7-ray pulsations in future Fermi data. 
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